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SUMMARY 


Theoretical analyses are presented of the effect of slip on the 
flow of a rarefied gas near a stagnation point and in a boundary layer 
on a flat plate. The results indicate that the stagnation pressure is 
increased because of the effect of slip but that there is a negligible 
effect on the flat-plate sk±n-friction coefficient in the range of 
application of the analysis. 


INTRODUCTION 


In the rarefied-gas-dynamic s regime it is known that the ordinary 
condition of no slip at a wall is altered and should be replaced by a 
relation between the wall shear and a slip velocity, usually taken to 
be given in the form indicated in equation (i|.) or equation ( 13 ) (refer- 
ence 1). It is also believed that the Navier-Stokes equations of gas 
dynamics are inadequate for accurate description of the flow field. 

At least three different systems of supposedly more adequate equations 
have been proposed in references 2 to 6 but their validity remains open 
to question. It is the purpose of this note to determine theoretically 
the rarefaction effect of slip in the boundary condition, using the 
NaviereStokes differential equations. As will be indicated below, it 
is probable that the rarefaction effect associated with the boundary 
condition produces terms of the same order as the still undetermined 
rarefaction effect associated with the differential equations. Hence, 
the present analysis probably furnishes only a partial improvement over 
the normal density treatment. The relatively sing^le boundary-layer case 
is considered first and then the more complicated stagnation-point flow. 

This work was conducted at the University of California under the 
sponsorship and with the fi n ancial assistance of the National Advisory 
Committee for Aeronautics and was under the supervision of Prof. R. G. 
Folsom, E. D. Kane, and S. A. Schaaf of the Department of Engineering. 
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SYMBCfLS 

b radius of curvature of nose of impact tube or radius of sphere 



f stagnation-point f\inction, defined in equations (16) and (I7) 

f^ stagnation-point function for no-slip case 

fj slip correction to stagnation-point function, equation (21;) 

K constant in equation (18) 

L molecular mean free path 

M Mach number 

p gas pressure 

p , mean hydrostatic stagnation pressure 

sx» 

p normal pressure in z-direction 

zz 

f p ) normal stagnation pressure 

\ st 

viscous correction to impact pressure, given by eqtiation (53) 
viscous correction to impact pressure for case of no slip 

r distance from z-axis 

Re Reynolds number 

Re^j Reynolds number based on x ^u^x^v) 

u tangential velocity component, that is, velocity component in 

x-direction; also radial velocity con^jonent, that is, 
velocity component in r-direction 
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Uq free-stream velocity past flat plate 

V normal velocity component, that is, velocity component in 

y-directionj also axial velocity component, that is, velocity 
con53onent in z-direction 

undisturbed free-stream velocity of gas 

X distance tangential to flat plate 

y distance normal to flat plate 

z distance normal to wall 

cc coefficient of sliding friction 

p velocity expansion parameter, given by equations (13) 

7 ratio of specific heats of gas 

6* boundary-layer displacement thickness 

6q* no-slip boundary-layer displacement thickness 

•q dimensionless independent variable, defined by equation {!$) 

X dimensionless parameter, defined by equation (20) 

[i. mechanical-viscosity coefficient 

V kinematic- viscosity coefficient ([x/p) 

i coefficient of slip, defined by equations (5) and (12) 

p gas density 

a Maxwell's reflection coefficient 

stream function, equations (2) 

\1/’ Blasius stream function 

'o 

slip correction to Blasius stream function 
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Subscripts: 

B condition at edge of boundary layer 

1 free-stream conditions 

o no-slip conditions 

BOUNDAHT-LAIER FLOW 


The stream function for incon5)ressible boundary layer on a flat 
plate is governed by the Blasiun differential equation (reference 7) 


a^T{f /.s 

3y ax ay ax ^ ^3 ^ ^ 

where 

* 

> ( 2 ) 

The corresponding boundary conditions are 

V = 0, y = 0 
u » Uq, y = to 

and for the raref ied-gas-dynami cs regime (reference l) 



By 


= 

dx 


t|.u, y-0 


(U) 


5 = 0.998 ^ ~ ^ L 
a 


( 5 ) 
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where a is the Maxwell reflection coefficient and L, the molecular 
mean free path. For the case of diffuse reflection | ~ L. In terms 
of the stream function \|c 


y = 0, y » 0 


L 


dy^ 



dy 


Uq, 7 = “ 


( 6 ) 


where equation ( 6 ) reduces to the ordinary density case if L = 0. Write 




Substitute into equations (l) and ( 6 ) and equate coefficients of L 
to 0 and obtain 




dy dx dy 3x 3^2 


= V 


ay- 


\lr = 0 and 
o 3y 


0 for y = 0 , and 




> (7) 


= Uq for y = « 


3x 3y 3y ^ 3y 9x 0 y .2 ^ 3x ^ 3 y .3 


= 0 and 


3ti 


> ( 8 ) 




^ „ for y = 0 , and = 0 for y = 

dy 3 y 2 3y 


Clearly the system of equation (?) is exactly the ordinary Blasius case 
and the solution is well-known. The solution of equation ( 8 ) may be 
directly verified to be given by 


^1 


5y 


(9) 
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so that \lr is immediately determined, to order L, from the known 
Blasius solution. The displacement thickness is given by 


6 * = 



^o ^ 


dy 



This may be written in the form 

^ - 0.867 (U) 


There is no change in the drag coefficient to the order of terms 
retained. This somewhat surprising result arises from the special 
nature of the assun^jtions made and it should not be inferred that there 
is no decrease in skin friction due to slip effects in a physical 
rarefied-gas flow. The mathematical expression indicated in equation (9) 
is a solution to the incompressible boundary-layer equation with a slip- 

flow boundary condition. Terms of the order ly^Re^^. as compared with 

terms retained have been systematically neglected in reducing the com- 
plete Navier-Stokes equations to the boundary-layer equation. The 
assumption of inconpressibility also limits the applicable range of 
Mach numbers. Thus the assumptions implicit in the use of equation (l) 
restrict the Mach and Reynolds numbers to roughly 


M < 0.5 
Re^. > IcA 


The effect of slip 
example, the term 

l/2 percent in the 
ever, is that this 
form 


on the boun dary-layer displacement thickness, for 
0.867My^^R^ in equation (U), is thus less than 


range of applicability. A more important fact, how- 
term is only one of several different terms of the 


Constant 




or 


Constant M^^Re^j. 
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which would presumably be obtained in a complete analysis, that is, a 
solution of the con 5 )lete Navier-Stokes equations for compressible flow. 
The present result with respect to skin friction thus seems to indicate 
only that the slip boundary condition, by itself, is initially of less 
importance than low Reynolds number modifications of boundary— layer 
theory in determining skin friction in a slightly rarefied gas. In this 
connection it should be noted that the somewhat cruder analyses of 
references 8 and ^ as well as the experimental data obtained at Berkeley 
and presented in references 10 and 11 (see fig. 1) all indicate a rela- 
tive decrease in skin friction due to slip effects at sufficiently low 
Reynolds numbers but only after the skin friction has first increased 
above the value it has at high Reynolds numbers. This is presumably 
associated with the breakdown of boundary-layer theory, that is, to the 
increasing importance of additional terms of the form Constant /^Re^. 

A further assun^ition implicit in the use of equation (l) is that 
rarefied-gas-dynamic terms, for example, Burnett terms, in the basic 
differential equations such as those discussed in references 2 to 6 may 
be neglected. Inspection of these equations suggests that additional 

terms of the form Constant /\^ would also arise in this connection. 
There are also additio n al terms in the boundary condition (reference U) 
but these do not affect the present results to the order of terms 
retained. In the stagnation-point-flow analysis presented in the next 
section, such effects have been neglectedj the complete Navier-Stokes 
equations have been solved to the order of terms retained, however. 


STAGNATION-POINT FLOW 


It is the purpose of this section to provide a theoretical basis 
for the interpretation of inqjact-tube measurements taken in a high- . 
velocity gas stream of sufficiently low density so that the effect of 
velocity slip on the reading of the instrument can no longer be ignored. 

Considered first is the problem of finding the flow in the vicinity 
of a forward stagnation point for a viscous fluid having the plane z = 0 
as a fixed wall. At the surface of the wall are prescribed the boundary 
conditions that the normal velocity is zero and that the tangential 
velocity - that is, the velocity of slip - is proportional to the 
velocity gradient in the normal direction. It is known from the kinetic 
theoiy of gases that the factor of proportionality | is essentially 
(reference l): 


I = 0.998 



L 


( 12 ) 
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where a is the Maxwell reflection coefficient and L, the molecular 
mean free path. At large distances from the wall the motion is essen- 
tially the same as that given by the potential flow. It follows that 
the boundary conditions for the problem are: 


For z = 0: 

u = I ^ and V = 0 
3z 


For z = “: 


(13) 


u = pr and v = -2p(z - 6*) ^ 


On assuming thAt the flow is steady, incompressible, and free from 
external forces, the forward stagnation-point-flow problem can be 
reduced to that of solving the third-order nonl i near ordinary differ- 


ential equation (reference 12), 

f I I I + 2ff • • - (f »)2 + 1 = 0 (Ik) 

satisfying certain conditions on the boundary with 

T) = 0/v)^/^z (15) 

u = prf'(il) (16) 

V = -2(pv)l/2f (ti) (17) 

and 

p = K - I p(p2r2 + v2) + Ij (18) 


The primes denote differentiation with respect to t). For the present 
problem, upon using equations (l5) to (17), the boundary conditions 
(equations (13)) become 

f(0) = 0 

f '(0) = Xf "(0) ^ 

(19) 


and 


f'(«) = 1 
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where 


X = l(p/v)^/^ = 0 . 998 ^-^ 

using equation (12). 



( 20 ) 


Now ii* fg denotes the function f for the nonslip case, one has 
(reference 12) 


fo"' + 2f^fo” - (fo’f ^ ^ ° 


( 21 ) 


and 


^o(O) = 0 
fo'(O) = 0 
fo’H = 1 



From Homann’s solution to f^ (reference 13) 

fo(“) = "H - 0.5576 
f^'»(0) = 1.317 } 

fo”(~) = 0 


( 22 ) 


(23) 


For \ small, one may write 

f = fp + Xf^ (2U) 

accurate to the first power of X. Then from equations (lij.), (19), (2l), 
(22), and (2i|) 



10 


MCA TN 2^68 


. and 


fl(0) = 0 

fl«(0) = fo"(0) V 


= 0 


( 26 ) 


also accurate to the first power of Differentiating equation (2l) 

with respect to t] once, one obtains 

f^iv 2f^^i.i = 0 (27) 

From equations (22) and (23) it is seen that f^ satisfies the conditions 


fo'(O) = 0 


( 28 ) 


If f]^ is replaced b 7 } equations (25) and (26) are reduced identi- 
cally to equations (2?) and (28), respectively. Therefore 

fq = fo' (29) 

is the solution to equations (25) and (26). Hence from equation {2k) 

f = f o + Xfo» (30) 

Thus, besides the parameter X, the function f of the slip case is 
expressed in terms of f^ and its derivative fQ* of the nonslip case. 
Making use of equation (30) and Homann's solution to fp (reference 13), 
one obtains in figure 2 the curves of f, f’, and f’* for X = 0.1. 

Eeferring to equation (50), X ~ 2 m/^ (see fig. 3). The corresponding 
curves for the nonslip case and the potential case are also shown in 
figure 2 for comparison. 

At the stagnation point r = z = v = 0 and p = Pg^. Substituting 
these into equation (l8) 
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p=t - - p(g),.o (31) 

where is the mean hydrostatic pressure, or the mean of the normal 

pressures across any three orthogonal planes, at the stagnation point. 

For ing)act-tube measurements the normal stress across .the wall p^g at 
the stagnation point is of special interest. For incon^jressible flow, 
in general (reference ll;), 

Pzz = -P + 2p (32) 

oz 

At the stagnation point 

^ Pzz ~ (Pzz)st 

From equations (32) and (3l) 

(Pzz)st = 

From equations (l5) and (1?) 

(3W 

From equation (30), upon using equations (23) and (22) 

f'(0) = 1.317>i (35) 

Substituting equation (3U) into equation (33), and using equation (35), 

(Pzz)st = "(K ^ 1'317 X 2ppx) (36) 

where the negative sign indicates compressive stress. 

The above equation is a solution to the Navier-Stokes equations with 
slip. So far no boundary-layer assumption has been made. 
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IMPACT TUBE 


Solution ( 36 ) will be applied to the case of an i mp act tube with a 
half-spherical nose, tested at subsonic speed at low pressure conditions, 
with its longitudinal axis placed parallel to the direction of the flow. 
Let the flow along the central streamline leading to the stagnation point 
decelerate isentropically to the outside seam of the boundaiy layer, 
which is presumed to form in the neighborhood of the stagnation point. 

It .is assumed that in the boundary layer in the neighborhood of the 
stagnation point the flow can be treated as viscous and incompressible 
and that outside the boundary layer the effects of viscosity can be 
neglected. 


At the point of the boundary layer where the stagnation streamline 

enters, that is, z = Zg, 

that from equation (I 8 ) 


r = 0 and f — ) = -26 (reference 1 $) so 

\BzJz=z^ 


K 


= ( 


% * ? 


+ 


(37) 


where pg and Vg are the pressure and velocity at the point Zg. The 

term in parentheses represents the i^^3act pressure p^^ recorded by an 

impact tube placed in a gas stream if the effect of compressibility in 
the boundary layer is zero. Since this was previously assumed for a 
subsonic flow, the impact pressure is given by the conventional isentropic- 
fldw equation (reference 16) 


Pk ■ Pb * 2 Pb V 




Substituting these into equation ( 36 ) 


(38) 



2pp(l + 1 . 3171 ) 


(39) 
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Tne constant p must be determined from, the velocity distribution 
at the seam of the boundary layer which surrounds the bocfy' at the stag- 
nation point. Since the flow outside the boundaiy layer is assumed 
inviscous and compressible, P can be deduced from the corresponding 
ideal subsonic compressible pattern. For a half-spherical impact tube 
of radius b, p has the form (reference 1?) 

^ = l(f - t V) (^ 0 ) 


where 


is the undisturbed free-stream velocity. 


Because of the presence of a boundary layer of average displacement 
thickness 6|*', the effective radius of curvature of the nose of the 
impact tube ^becomes b + 6*. For viscous incompressible flow passing a 
sphere, the ^average displacement thickness of the boundaiy layer 

of the nonsii.p case can be expressed as (reference 13) 


where 


Similarly 



V = "H - fo(“) 


pV.b 

Re = 

V- 


6 * _ T]* 



m 


ih2) 

ih3) 


with 


-n* = T) - f (~) 


From equation (30), upon using equations (23) and (22), 
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so that 


Ti* = 0.5576 - \ 


and 


^ ^ (0.1i55 - 0.8161) 

Therefore equation (UO) becomes 

(29 31 j.2\ 

^ 2b L , 0.1i55 - 0.8l6\ 

\ 


(115) 

(116) 


(li7) 


Substituting this into equation (39) 



I L 7 - 1 „ 2\y-^ 


(f ~ i 


o.li55 - o.8i6\\ 


0i8) 


j/ie 


Consider next the expression for X. Substituting equation (li8) 
into equation (20) and for the present, since X is small, neglecting 
the slip correction to the effective radius of the iinpact tube. 


\ cr JU + 0.1i55/\|^y ^ 

Upon using the relation (reference 18) 

L/h = 1.255f7 M/Re 

with y = I.I 1 O 5 and a = 1, equation (h9) becomes 


(li9) 


( 50 ) 
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X 


'i - 

+ 0.li55/#i/ 



(51) 


The vaxiation of X with the rarefaction parameter M^/Re for various 
values of the Mach number M-j_ is shown in figure 3» Substituting 
equation (5l) into equation (ij.8) 



If denotes the viscous correction to the impact-tube pressure, 
from equation (I4.8) 


P.. = 


p,Vi 


(f -f 


1 + 


0.1;5 5 - 0 .816X 

\f^ 


(53) 


For the nonslip case X = 0^ Pp, - Ppo* From equation (53) 



(5it) 


so that, upon using equation (5l) 
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(l - 0.252M^2^^/2 ji 

1 + 2.63 ) ViV9 

Pg (l + O.li^^AlRi) ^ ^ 

fl - 0 . 252 Mn 2 )^/^ II, 




For Mj^Re .gmall and Mj small the above expression is larger than 

■unity. It follows that the effect of slip increases the viscous cor- 
rection to the inpact-tube measurements. The curves of 

against Kj^^Re, l/Re, and Re for various values of are given 

in figures I 4 ., and 6, respectively. 


Define 



(56) 


Upon using equations (56), (53), (52), and (5l), the -variations of Cjj_ 

and Cm with Re and l/Re for Mq = O.i; and 0.8 are shown in fig- 
^*"0 

ures 7 and 8 , respectively. 


SLIP FLOW PAST A SPHERE 


It is interesting to note that at the stagnation point the nor mal 
pressure across the wall increases with the coefficient of slip, 
although at the same point the mean hydrostatic pressure, that is, the 
mean of the normal press'ures across any three orthogonal planes, 
decreases. This is due to the fact that for the slip case the normal 
velocity gradient (dv/dz) at the stagnation point is no longer zero. 
It follows that the effect of slip increases the viscous correction 
(reference l5) to the impact-tube measurements. A similar result is 
obtained for 'the case of the slow translation of a sphere through a. 
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viscous fluid vrith slip from the well-known solution given hy Basset 
(reference 19 or reference li;, pp. 603 - 60 J 4 .) . In his case the inertia 
terms are neglected. 

For the case of the slow translation of a ^here through a viscous 
fluid with slip the i mpact pressure p^ at the forward stagnation point 

has been obtained by Basset (reference 19) j who neglected the inertia 
terms which are small in comparison with the viscous and pressure terms. 
If b is the radius of the sphere, the velocity, a, the coef- 

ficient of sliding friction, and p , the hydrostatic pressure, then 


p .p 

^st ^1 2 b T + ^ 


On using the relation (reference I 8 ) 


ab 


i - i.253Vr(^)s| 

equation (5?) becomes, for y ~ l.ii 05 and a = 1, 


(57) 


(58) 



+ 


3 1 + 8.88M]^/Re 

2 b 1 + 


(59) 


The viscous correction to the impact pressure is 

3 pVi 1 + 8.88M3^/Re 

" 2 ~ 1 + li.lUiM^Re 


For the nonslip case a = 00 ^ 


from equation (57) 


Therefore 



3 ^"^1 
2 b 


^ _ 1 + 8.88M]^/Re 
1 + ii.ijl4lI]yRe 


(60) 


( 61 ) 


( 62 ) 
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As a comparison the curve of "tiie Basset sphere against 

M^^^Re is shovm in figure ij.. 


CONCLUSIONS 


A method has been developed for determining some of the rarefaction 
effects in slip-regime flow for the flat plate, a wall., and a hemi- 
spherical impact tube, using a reduced form of the Navie3>-Stokes equa- 
tions and a slip boundary condition. With this method it was found 
that: 


1. The effect of slip on the viscous flow about an inqpact tube is 
to increase the viscous correction to the ideal imp act pressure. 

2. The effect of slip on the flat plate is to reduce the boundary- 
layer displacement thickness. There is no effect on the drag coef- 
ficient to the order of terms retained in the analysis. 


University of California 

Berkeley, Calif., October 20, 1950 
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Figure 2.- Functions of forward stagnation-point flow, u ° prf'(T)). 
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Figure 7.- Viscous 
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